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ABSTRACT 

We present nearly simultaneous Chandra and NuSTAR observations of two actively star-forming galax¬ 
ies within 50 Mpc: NGC 3256 and NGC 3310. Both galaxies are significantly detected by both Chandra 
and NuSTAR, which together provide the first-ever spectra of these two galaxies spanning 0.3-30 keV. The 
X-ray emission from both galaxies is spatially resolved by Chandra ; we find that hot gas dominates the 
E < 1-3 keV emission while ultraluminous X-ray sources (ULXs) provide majority contributions to the emis¬ 
sion at E > 1-3 keV. The NuSTAR galaxy-wide spectra of both galaxies follow steep power-law distributions 
with T « 2.6 at E > 5-7 keV. Using new and archival Chandra data, we search for signatures of heavily ob¬ 
scured or low luminosity AGN. We find that both NGC 3256 and NGC 3310 have X-ray detected sources 
coincident with nuclear regions; however, the steep NuSTAR spectra of both galaxies restricts these sources to 
be either low luminosity AGN (L 2-10 kev/^Edd ~ 10“ 5 ) or non-AGN in nature (e.g., ULXs or crowded X-ray 
sources that reach L 2-10 keV ~ 10 4(1 erg s _1 cannot be ruled out). Combining our constraints on the 0.3-30 keV 
spectra of NGC 3256 and NGC 3310 with equivalent measurements for nearby star-forming galaxies M83 and 
NGC 253, we analyze the star-formation rate (SFR) normalized spectra of these starburst galaxies. The spectra 
of all four galaxies show sharply declining power-law slopes at energies above 3-6 keV primarily due to ULX 
populations. Our observations therefore constrain the average spectral shape of galaxy-wide populations of 
luminous accreting binaries (i.e., ULXs). Interestingly, despite a completely different galaxy sample selection, 
emphasizing here a range of SFRs and stellar masses, these properties are similar to those of super-Eddington 
accreting ULXs that have been studied individually in a targeted NuSTAR ULX program. We also find that 
NGC 3310 exhibits a factor of ft3—10 elevation of X-ray emission over the other star-forming galaxies due 
to a corresponding overabundance of ULXs. We argue that the excess of ULXs in NGC 3310 is most likely 
explained by the relatively low metallicity of the young stellar population in this galaxy, a property that is 
expected to produce an excess of luminous X-ray binaries for a given SFR. 

Subject headings: galaxies: individual (NGC 3256 and NGC 3310) — galaxies: active — galaxies: starburst 
— galaxies: star formation — X-rays: galaxies 
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With the launch of the NuSTAR X-ray observatory (Harri¬ 
son etal. 2013), we now have a first-ever means for study¬ 
ing populations of normal galaxies (not hosting luminous ac¬ 
tive galactic nuclei [AGN]) in the ft 10-30 keV bandpass, 
an energy regime which is expected to be dominated by 
bright X-ray binaries and low-level contributions from diffuse 
inverse-Compton emission (e.g., Persic & Rephaeli 2002; 
Wik etal. 2014a). The spectra of X-ray binaries at these en¬ 
ergies betray unique information about their accretion states 
as variations in accretion disk, reflection, and coronal compo¬ 
nents undergo strong spectral pivots above ftlO keV (see, e.g.. 
Done etal. 2007), a regime that was not previously probed 
by Chandra and XMM-Newton for sources outside the Local 
Group. Prior to NuSTAR, studies of the hard X-ray spectra of 
X-ray binaries was limited to objects in the Milky Way (MW) 
and Magellanic Clouds (e.g., via RXTE\ see Remillard & Mc- 
Clintock 2006 for a review), providing limited information 
about how the accretion states of these populations vary with 
galaxy properties (e.g., the starburst environment). 

With Chandra and XMM-Newton, the lower-energy 
<10keV emission from local galaxies has been studied in 
detail for a variety of galaxy types (see, e.g., Fabbiano et al. 
2006 for a review), and ultradeep Chandra observations have 
extended these studies to very distant galaxy populations with 
look-back times that span the vast majority of cosmic history 
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Fig. 1.— Physical properties of NuSTAR starburst galaxies sample. ( a ) SFR versus M* for the full sample with regions of constant specific SFR 
(sSFR = SFR IM+) highlighted. ( b ) sSFR versus distance for the full sample. For comparison, the location of the Milky Way (MW) has been indicated. 
Our sample consists of galaxies with sSFRs that range from being comparable to those of the MW to a factor of «20 times higher. We therefore expect that the 
X-ray emitting point-source populations in our sample will be skewed towards including richer populations of HMXBs compared to the MW (e.g., Lehmer etal. 
2010; Mineo etal. 2012a). NGC 3256 and NGC 3310 are amongst the highest sSFR galaxies in our sample. 


(e.g., Basu-Zych etal. 2013a). We are detecting X-ray emis¬ 
sion from galaxy populations out to z « 4 where the observed- 
frame Chandra 0.5-8 keV band corresponds to rest-frame en¬ 
ergies of 2.5-40 keV, the most sensitive regimes of the NuS¬ 
TAR bandpass for extragalactic binary work. Results from 
such studies have prompted the development of accreting bi¬ 
nary population synthesis models, which have successfully 
been employed to construct a self-consistent framework that 
describes the observed evolution of X-ray binary populations 
from cosmic dawn to today (e.g., Fragos etal. 2013a,b). Key 
to these observational interpretations and model predictions, 
however, is knowledge of the X-ray spectra at x) 10 keV for 
galaxy populations that cover a broad range of physical prop¬ 
erties (e.g., Fragos etal. 2013a; Kaaret etal. 2014). There¬ 
fore, a next step to improving our understanding of X-ray 
binary populations is to constrain the distributions and duty 
cycles of their accretion states and measure how their result¬ 
ing xl 10 keV spectral contributions vary with galaxy physi¬ 
cal properties (e.g., star-formation rate [SFR] and stellar mass 
[M*]). 

We are conducting a joint NuSTAR and Chandra survey 
of six far-infrared bright starburst galaxies (NGC 253, M82, 
M83, Arp 299, NGC 3256, and NGC 3310) with the goal 
of quantifying the dominant processes that contribute to the 
X. 10 keV emission. Figure 1 displays the SFR versus M* and 
specific SFR (sSFR = SFR/M*) versus distance planes for the 
full galaxy sample. Our key goals are to: (1) identify the ac¬ 
cretion states of the most luminous X-ray binary populations 
in these galaxies and characterize the «0.3-30 keV spectral 
energy distributions (SEDs) as a function of the galaxy phys¬ 
ical properties (e.g., SFR and M*); (2) search for heavily ob¬ 
scured AGN that may be present in these actively star-forming 
galaxies; and (3) constrain the nature of inverse Compton 
emission associated with particle accelerations in starburst 
flows. 

As a pilot program, we studied with NuSTAR and Chan¬ 
dra the nuclear region and galaxy-wide X-ray emission of the 


nearby starburst galaxy NGC 253 (Fehmer etal. 2013; Wik 
etal. 2014a). These studies showed that there was no evi¬ 
dence for powerful buried AGN activity in the nucleus, and 
that the galaxy-wide X-ray emission above 10 keV is domi¬ 
nated by a few ultraluminous X-ray sources (UFXs) with mi¬ 
nority contributions from lower luminosity X-ray binaries that 
have NuSTAR colors similar to Galactic black hole binaries in 
intermediate accretion states. The galaxy-wide X-ray spec¬ 
trum of NGC 253 steepens at energies above «6 keV, consis¬ 
tent with the spectra of other UFXs studied by NuSTAR (see, 
e.g., Bachetti etal. 2013; Rana etal. 2014; Walton et al. 2013, 
2014), signaling a possible dominance of super-Eddington ac¬ 
creting objects in the overall starburst galaxy spectra. 

More recently, we have executed joint NuSTAR and Chan- 
dra/XMM-Newton observations of M83 (Yukita etal. in¬ 
preparation) and Arp 299 (Ptak et al. 2014), galaxies with sS¬ 
FRs that are comparable to and «20 times higher than that of 
the MW, respectively. The brightest X-ray sources in M83 
have NuSTAR colors similar to those of NGC 253, albeit 
with fewer UFXs. Arp 299 has x( 10 keV emission domi¬ 
nated by a single Compton-thick AGN that outshines the other 
X-ray emitting components of the galaxy. Given the relatively 
large distance (rs 40 Mpc) to Arp 299, and its correspondingly 
smaller angular extent, it was not possible to spatially resolve 
the X-ray emitting components with NuSTAR. Therefore, the 
spectral properties of the non-AGN components were poorly 
constrained at energies above «10 keV, leaving us with little 
knowledge of the high-energy spectral shape for high sSFR 
galaxies (see Fig. lb). 

In this paper, we continue to investigate the 0.3-30 keV 
SEDs of starburst galaxies by studying two powerful starburst 
galaxies in the local universe: NGC 3256 and NGC 3310. 
The sSFRs of these galaxies are the second and third high¬ 
est of our sample (see Fig. 1), and are only exceeded by 
Arp 299. Both NGC 3256 and NGC 3310 have been studied 
extensively across the full electromagnetic spectrum, and nei¬ 
ther system exhibits compelling evidence for harboring pow- 




NuSTAR and Chandra View of NGC 3256 and NGC 3310 


3 


erful AGN (see §2 for a discussion). The sSFRs of these 
galaxies are «0.4-0.8 Gyr 1 , a range expected to host X-ray 
binary populations that are dominated by high-mass X-ray 
binaries (HMXBs; e.g., Colbert etal. 2004; Lehmer etal. 
2010). Such high-sSFR galaxies are representative of the typ¬ 
ical galaxy populations on the “main sequence” at z ~ 1-2 
(e.g., Karim etal. 2011). As such, this study will provide 
SED constraints on high-sSFR galaxies that can be used to 
inform X-ray studies of high-redshift galaxy populations at 
these redshifts, which require informed K-corrections (e.g., 
Lehmer etal. 2008; Basu-Zych et al. 2013a; Kaaret 2014). 

Throughout this paper, we assume distances of 
D = 35.4 Mpc and 19.8 Mpc (Sanders et al. 2003) and column 
densities of Ah = 9.6 x 10 20 cm -2 and 1.1 x 10 20 cm -2 
for NGC 3256 and NGC 3310, respectively (Dickey & 
Lockman 1990). Star-formation rates and stellar masses 
quoted throughout this paper were calculated assuming a 
Kroupa (2001) initial mass function. Star-formation rates 
were calculated using equation 6.11 of Calzetti (2013), 
making use of UV and mid-IR data from GALEX and Spitzer, 
respectively. Stellar masses were estimated following the 
prescription outlined in Appendix 2 of Bell etal. (2003), 
making use of B — V colors from RC3 (de Vaucouleurs 
etal. 1991) and K-band luminosities based on 2MASS 
photometry (Jarret etal. 2003) and the adopted distances. 
These assumptions imply SFR ss 36 M 0 yr 1 and 6 Mq yr 1 
and M*«9x 10 10 Mq and 9 x 10 9 M 0 for NGC 3256 
and NGC 3310, respectively. Quoted errors associated with 
spectral fits represent 90% confidence intervals. 

2. THE GALAXIES 
2.1. NGC 3256 

NGC 3256 (35.4 Mpc) is a major-merger system composed 
of two gas-rich galaxies with comparable masses (total M* r; 
10 11 M 0 ) that are in a nearly coalescent phase of the merger 
(Lfpari etal. 2000; Alonso-Herrero etal. 2002). Two discrete 
nuclei have been identified «5 arcsec apart (850 pc; Zenner 
& Lenzen 1993; Norris & Forbes 1995) and remarkably long 
(200 kpc) H I tidal features have been identified (e.g., Gra¬ 
ham etal. 1984; English etal. 2003). The system has the 
highest 8-1000 pm infrared luminosity (Lir « 4x 10 11 L 0 ; 
Sanders etal. 2003) for galaxies that reside within z A 0.01, 
and studies at radio (e.g., Norris & Forbes 1995), submm 
(e.g., Sakamoto etal. 2006, 2014), infrared (Graham etal. 
1984; Doyon etal. 1994a, 1994b; Moorwood & Olivia 1994; 
Lira etal. 2008; Alonso-Herrero etal. 2013), optical (e.g., 
Heckman et al. 2000; Lfpari et al. 2000; Alonso-Herrero et al. 
2002), and UV (e.g., Kinney etal. 1993; Leithereretal. 2013) 
wavelengths all support a scenario in which the majority of 
the galaxy’s power output is produced by starburst-related ac¬ 
tivity. 

X-ray studies with Chandra and XMM-Newton have shown 
that the 0.3-10 keV spectrum of NGC 3256 can be modeled 
well as a multiphase thermal plasma with harder power-law 
emission from a population of discrete point sources, the ma¬ 
jority of which are ULXs (see, e.g.. Lira etal. 2002; Jenkins 
etal. 2004; Pereira-Santaella etal. 2011). Although both nu¬ 
clear regions have been detected by Chandra with L 2-10 keV ~ 
10 4() erg s -1 , the nature of the X-ray emission is unclear, as no 
obvious AGN signatures have been found (Lira etal. 2002). 
Detailed assessment of the kinematics of the molecular gas, 
via ALMA and SMA observations, have concluded that both 
nuclei contain powerful outflows, but with differing character- 
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Fig. 2.— Relative NuSTAR (black lines) and Chandra {gray lines) observa¬ 
tional coverage for NGC 3256 {top) and NGC 3310 {bottom). For clarity, we 
have annotated the stalling date of the observational epoch. The apparently 
broken up NuSTAR observational intervals are due primarily to Earth occul- 
tations and passages through the SAA, which result in observing efficiencies 
of s=:51-58% on average. 

istics (e.g., Cicone etal. 2014; Sakamoto etal. 2014; Emonts 
etal. 2014). The northern nuclear outflow is of wide breadth 
and has a mean velocity of a few 100 km s -1 and is likely 
driven by a starburst superwind. The southern outflow is well- 
collimated and faster ( ^ 1000 km s -1 ) and is plausibly driven 
by an AGN that has recently gone dormant (e.g., Moran et al. 
1999; Neff et al. 2003; Alonso-Herrero et al. 2012). 

2.2. NGC 3310 

NGC 3310 is a closer (19.8 Mpc) starburst system with 
a SLR that is a factor of ~10 times lower than NGC 3256 
(Lir ss 4 x 10 10 L 0 ; Sanders etal. 2003); however, given 
its relatively low stellar mass (M* 10 10 M 0 ), this system 

has a comparable and perhaps larger sSLR than NGC 3256 
(see Fig. 1). The system contains evidence for a minor- 
merger that has triggered a vigorous, young («2.5-5 Myr 
old) circumnuclear star-forming ring with a ~20-30 arcsec 
(r;2 kpc) diameter (e.g., Balick & Heckman 1981; Dfaz etal. 
2000; Elmegreen etal. 2002; Chandar etal. 2005; Miralles- 
Caballero et al. 2014). The most recent minor merger was in¬ 
ferred to occur in the last r;30 Myr (Elmegreen etal. 2002; 
de Grijs etal. 2003a), and is likely to be one of several 
past interactions with small metal-poor dwarf galaxies (e.g., 
Wehner et al. 2006), as evidence for starburst activity as old as 
100 Myr is apparent in the star clusters (Meurer 2000). Such a 
series of minor-mergers could potentially provide a significant 
increase to the mass build-up of the bulge (Miralles-Caballero 
et al. 2014) and is likely to have modified the gas-phase abun¬ 
dances, as there is evidence that the circumnuclear metallicity 
is lower than is typical for spiral galaxies with morphologies 
similar to NGC 3310 (e.g., Pastoriza et al. 1993). 

Similar to NGC 3256, the Chandra and XMM-Newton spec¬ 
tra of NGC 3310 can be modeled well by a multiphase hot gas 
plus a power-law component; however, the plasma tempera¬ 
ture is cooler and the power-law slope is harder than those 
found for NGC 3256 (Jenkins etal. 2004). The hard power- 
law component is dominated by a population of «9 ULXs that 
are distributed within the circumnuclear star-forming ring and 
along the spiral arms. HST STIS observations of the nuclear 
region have uncovered a bright rT).5 arcsec (r;40 pc) diameter 
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Fig. 3.— Three-color adaptively smoothed Chandra images of NGC 3256 and NGC 3310 with 4-25 keV NuSTAR intensity contours overlaid ( dashed curves). 
The three bandpasses, 0.3-1 keV (red), 1-2 keV (green), and 2-7 keV (blue) were chosen to highlight the diffuse components from hot gas and unresolved 
point sources, as well as absorption in the bright point-source populations. A 10 arcsec bar has been added to each figure for reference. The annotated source 
numbers, ordered by increasing right ascention, highlight bright point sources that are studied in §4.3 and presented in Table 2. The locations of galactic nuclei 
are indicated with magenta crosses. For NGC 3256, the northern (N) and southern (S) nuclei are labeled. With Chandra, we detect the northern nucleus of 
NGC 3256 (source 4, left) and the nuclear region of NGC 3310 (source 3, right). The NuSTAR FWHM is «18 arcsec, which is sufficient for resolving some 
structure in the nearer NGC 3310; however, the NuSTAR emission for NGC 3256 is consistent with a single source PSF. 


central star-forming region with kinematic properties consis¬ 
tent with a circularly rotating disk (Pastorini et al. 2007). The 
nuclear region is detected by Chandra as a point-like sourceFl 
and the luminosity of this region, L 2-10 keV ~ 10 4<) erg s _1 , is 
within the range of other ULXs in the galaxy. The Chan¬ 
dra spectrum shows that the nuclear X-ray source has a flat 
power-law spectral slope, with some evidence (at the 2a con¬ 
fidence level) for an Fe Ka line—two properties that indicate 
there may be a hidden (obscured) AGN in the galaxy nucleus 
(Tzanavaris etal. 2007). 

3. OBSERVATIONS AND DATA REDUCTION 

In this investigation we make use of nearly simultaneous 
observations with NuSTAR and Chandra to explore the full 
0.3-30 keV bandpass properties. In later sections, we make 
use of archival Chandra data to support our analyses; how¬ 
ever, we focus the current investigation on the new nearly 
simultaneous observations. We obtained our NuSTAR and 
Chandra observations of NGC 3256 and NGC 3310 over sin¬ 
gle epochs beginning on 2014 Aug 24 and 2014 Jun 11, re¬ 
spectively. Figure 2 shows the relative observing schedules 
for the NuSTAR and Chandra exposures. Due to Earth oc- 
cultations and passages through the South Atlantic Anomaly 
(SAA), the NuSTAR on-target observations were carried out 
at «51-58% efficiency over «3-4 day periods. The cumula¬ 
tive NuSTAR good-time-interval exposures were 184 ks (Ob- 
sID 50002042) and 141 ks (ObsID 50002040) for NGC 3256 
and NGC 3310, respectively. The Chandra observations were 
continuous and lasted 16 ks (ObsID 16026) and 10 ks (ObsID 
16025) for NGC 3256 and NGC 3310, respectively. 

14 The nuclear region is also detected by XMM-Newton\ however, it is 
clear that confusion with four other Chandra- resolved ULXs of comparable 
brightness makes it difficult to uniquely interpret those data. 


The NuSTAR data were reduced using HEASoft v6.15, 
NuSTAR Data Analysis Software (NuSTARDAS) vl.7, and 
CALDB version 20130320. We processed level 1 data to level 
2 products by running nupipeline, which performs a va¬ 
riety of data reduction steps, including (1) filtering out bad 
pixels, (2) screening for cosmic rays and observational inter¬ 
vals when the background was too high (e.g., during passes 
through the SAA), and (3) projecting accurately the events to 
sky coordinates by determining the optical axis position and 
correcting for the dynamic relative offset of the optics bench 
to the focal plane bench due to motions of the 10 m mast that 
connects the two benches. 

Due to the angular extents of NGC 3256 and NGC 3310 
(major-axes 2 a = 3.8 arcmin and 3.1 arcmin, respectively), 
the Chandra exposures for both galaxies were conducted us¬ 
ing single ACIS-S pointings with the approximate position 
of the galactic centers set as the aimpoints. For our data re¬ 
duction, we used CIAO v. 4.6 with CALDB v. 4.6.1.1. We 
reprocessed our events lists, bringing level 1 to level 2 using 
the script chandra_repro, which identifies and removes 
events from bad pixels and columns, and filters events lists to 
include only good time intervals without significant flares and 
non-cosmic ray events corresponding to the standard ASCA 
grade set (grades 0, 2, 3,4, 6). We constructed an initial Chan¬ 
dra source catalog by searching the 0.5-7 keV images with 
wavdetect (run with a point spread function [PSF] map 
created using mkpsfmap), which was set at a false-positive 
probability threshold of 2 x 10~ 5 and run over seven scales 
from 1-8 (spaced out by factors of \J~2 in wavelet scale: 1, 
v/2, 2, 2^/2, 4, 4-\/2, and 8). Using these initial Chandra 
source catalogs, we constructed source-free light curves and 
searched for flaring intervals that were 5er above the back¬ 
ground level. No such intervals were found, and our observa- 




NuSTAR and Chandra View of NGC 3256 and NGC 3310 


5 




0.5 1 2 5 10 20 

Energy (keV) 


Fig. 4.— Joint nearly simultaneous Chandra (gray) and NuSTAR (red for FPMA and blue for FPMB) spectra for the full galactic extents of NGC 3256 (left) 
and NGC 3310 (right). The diffuse emission that is spatially isolated with Chandra was used to constrain two-temperature thermal (apec) models for both 
galaxies. For NGC 3256, an additional obscured thermal component («1 keV) was added to account for hot gas emission associated with the nuclear starburst. 
Absorbed single and broken power-law models that account for the Chandra- detected binaries were both tested for each galaxy, and the resulting residuals are 
displayed, and goodness-of-fit values are annotated, in the bottom two panels of each plot. We find that a single power-law is sufficient to fit the non-thermal 
emission from NGC 3256, while a broken power-law is preferred in the case of NGC 3310. We also display the contributions from nuclear components of 
these galaxies, with extrapolations to the NuSTAR band. For NGC 3256, this includes contributions from two nuclear regions associated with two galaxies in the 
merger: the northern nucleus contribution (green curve) and the estimated upper limit to the southern nucleus contribution (orange curve). The southern nucleus 
was not formally detected in our 2014 observation, and the upper limit spectrum here indicates the average Chandra spectrum from deeper observations in 2000 
and 2003 scaled to the upper limit flux of our 2014 observation. We note that the southern nuclear region of NGC 3256 and the nuclear region of NGC 3310 both 
have flat X-ray spectra and are candidate low luminosity AGN (see §4.3 for details). 


tions were deemed to be sufficiently cleaned. 

In Figure 3, we show the adaptively smoothed, false color 
Chandra images of NGC 3256 and NGC 3310, with NuS¬ 
TAR 4-25 keV intensity contours overlaid. It is apparent from 
the Chandra data that diffuse emission at E < 1 keV, poten¬ 
tially due to hot gas and unresolved point sources, and multi¬ 
ple bright point sources at all energies provide dominant con¬ 
tributions to the Chandra- detected emission for both galax¬ 
ies. The NuSTAR PSF has an 18 arcsec full width half max¬ 
imum (FWHM) core with a 58 arcsec half-power diameter 
(Harrison etal. 2013). Given the X-ray extents of the galax¬ 
ies (60-70 arcsec), NuSTAR is unable to strictly resolve the 
point-source populations that dominate the E > 3 keV emis¬ 
sion probed by NuSTAR ; however, some evidence for spatial 
extent following the point-source distribution is apparent in 
NGC 3310 (see §4.3 below). 

4. ANALYSIS AND RESULTS 

As discussed in §1, our primary goals are to measure 
and assess the contributing components to the broad-band 
0.3-30 keV SEDs of the two powerful starburst galaxies 
NGC 3256 and NGC 3310, and constrain the underlying AGN 
activity in these galaxies. To achieve these goals, we first 
characterized the 0.3-30 keV spectra across the full extents 
of both galaxies. We then used the high spatial resolution of 
Chandra to measure the 8 keV contributions from bright 
resolved point-sources and diffuse emission from hot gas and 
unresolved point-sources. Finally, using sensible extrapola¬ 
tions of the Chandra component spectra to £ 8 keV, we infer 
the relative contributions AGN candidates would make to the 


NuSTAR galaxy-wide spectra. In the sections below, we de¬ 
scribe each of these procedures in detail. 

4.1. Galaxy-Wide Spectral Analyses 

We began by extracting the nearly simultaneous Chandra 
and NuSTAR spectra over the full extents of both galaxies. 
For NGC 3256 and NGC 3310, we extracted on-source spec¬ 
tra using circular apertures with radii r src equal to 80 arcsec 
and 90 arcsec, respectively. These apertures were chosen to 
encompass the entire optical extents of the galaxies and reach 
20 arcsec (roughly the NuSTAR FWHM) beyond the most off¬ 
set point-sources that were clearly detected with Chandra. 
For each galaxy, background spectra were extracted using 1-5 
circular apertures located in source-free regions. For Chan¬ 
dra, the background extraction regions had radii «50-70 arc¬ 
sec and were chosen by eye to be placed in a pattern surround¬ 
ing the on-source extraction region. For NuSTAR, background 
regions were chosen more carefully to properly account for 
the spatial background gradients that arise primarily from the 
“aperture” background component, which contains cosmic 
X-ray background stray light that shines directly onto the de¬ 
tectors from a RilM-deg annular region (see, e.g.. Appendix A 
of Wik etal. 2014b). Following the procedure described in 
Appendix B of Wik et al. (2014b), we constructed background 
maps for the FPMA and FPMB modules in the 3-20 keV 
band, an energy regime sensitive to variations in the aperture 
background component. For each background map, we iden¬ 
tified 1-4 regions with sizes equal to the on-source aperture 
that had background levels comparable to those expected to be 
present within the source extraction regions. The closest edge 
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TABLE 1 

Best fit parameters for full-galaxy spectra of NGC 3256 and NGC 3310. 


Parameter 

Unit 

NGC 3256 

NGC 3310 

Chandra only fit:* point-source subtracted diffuse emission 

Ah,G al ... 


10 22 cnT 2 

0.096 

0.011 

Wi. 

Ah ,2 ■ ■ ■ ■ 

kT 2 . 


keV 

10 22 cm -2 
keV 

r\ 70 +O.O 6 
U-JJ -0.04 

0 74+0.09 

U - /4 -0.04 

0 on+0.06 
u - vu -0.08 

0.2 ± 0.1 

06 X):i 

°- 28 $ol 


Chandra + NuSTAR galaxy-wide fit:^ diffuse emission 

plus non-thermal point 

sources 

Ah ,3 — 


10 22 cnT 2 

115!‘ 28 


m . 

Ah.pl ■ ■ ■ 

Tj . 


keV 

10 22 cnT 2 

1 o +2 - 4 

3.0+JI 

n oq+0.30 
u< ”-0.56 

2.07 4 fJ';„ 

^break • ■ ■ 


keV 


4 70+ 0 - 72 

7 -1.21 

r 2 . 



2.64 ±0.20 
(1.77 ± 0.20) x 10~ 4 
(1.29 ^’jq) x 10“ 3 
4.0^ 6 x 10~ 2 

2 65 4-010 
z ’ dj -0.14 

(1.15 ±0.32) x 10 -4 
2 -85Lg;|| x IffA 

AkT 2 ■ ■ • ■ 
A k T 3 .... 



Apl . 


photons keV -1 cm -2 s -1 at 1 keV 

4.6+j | x 10^ 

7.4f 2 | X 10^ 

X?ed •••• 



1.09 

0.98 

V . 


(degrees of freedom) 

159 

192 

/x (10-30 keV) .. 

10“ 13 ergs cm -2 s -1 

1.33±0.03 

5.9±0.4 

L x (10-30 keV) ... 

10 40 ergs s -1 

1.99±0.05 

2.8±0.2 


Note—All quoted errors are at the 90% confidence level. 

‘Fits to the point-source subtracted diffuse Chandra spectra were performed using XSPEC model tbabSQai* (apeci+tbabS2*apeC2+tbabspL*pow) to 
obtain values of T \, 7), and Afi 2 . In this process the power-law had a fixed photon index of T = 1.8 and all thermal models assumed solar abundances. Values of 
T\, 73. and /\'n .2 were subsequently fixed when fitting the Chandra + NuSTAR galaxy-wide spectra, where their normalizations were free to vary, 
t The galaxy-wide Chandra + NuSTAR spectra were fit using XSPEC models tbabSGal* (ape c 1 +1 b a b s 0 * a p e c 2+1 b a b s 3 * a p e c 3+1 b a b .0 pl *bknp c ) and 
tbabSQai* (apeci+tbabS2*apeC2+tbabspL*pow) for NGC 3256 and NGC 3310, respectively. 


of each background region was located > 1.2 x r src away from 
the central positions of each galaxy to avoid contamination 
from the galaxy itself. After choosing appropriate regions, 
the Chandra and NuSTAR on-source and background spectra 
were extracted using the specextract and nuproducts 
tools, respectively. 

In Figure 4, we show the 0.3-30 keV galaxy-wide spec¬ 
tra for both galaxies. As discussed in §3, from the Chan¬ 
dra data (see Fig. 3), it is clear that the X-ray emission from 
both galaxies can be broadly characterized as consisting of 
diffuse emission (from hot gas and unresolved point sources) 
and bright point sources (e.g.. X-ray binaries and possi¬ 
bly AGN). Previous investigations with Chandra and XMM- 
Newton have found that the diffuse components of these 
galaxies can be modeled well using two-temperature thermal 
spectra, with NGC 3310 having cooler temperatures ( kT = 0.3 
and 0.6 keV) than the more powerful NGC 3256 (kT = 0.6 and 
0.9 keV), which may also contain an additional absorbed hot 
(«3.9 keV) component (e.g.. Lira etal. 2002; Jenkins etal. 
2004). The brightest point sources in these galaxies are ULXs 
with nine sources spanning Lx ~ (2-10) xlO 39 erg s' 1 for 
NGC 3310 (Smith etal. 2012) and roughly a dozen sources 
with L x ~ (2-60) xlO 39 erg s” 1 in NGC 3256 (Lira etal. 
2002). These ULXs dominate the galaxy-wide emission from 
«2-8 keV and are expected to provide majority contributions 
to the non-AGN emission in the ~8-30 keV bandpass. We 
therefore chose to model the broad-band 0.3-30 keV spectra 
of the galaxies using the sum of thermal and power-law com¬ 
ponents to account for the hot gas and ULXs, respectively. 

For each galaxy, we tested for the presence of 2-3 ther¬ 
mal components, in which the absorption of each compo¬ 
nent increases with increasing temperature. Such a trend has 
been seen in spatially resolved starburst galaxies like M82, 


NGC 253, and past studies of NGC 3256, in which the hottest, 
most obscured components are located in the nuclear regions 
and the cooler, unobscured plasmas extend to larger galaxy¬ 
wide scales (e.g., Strickland etal. 2000; Pietsch etal. 2001; 
Lira et al. 2002; Strickland & Heckman 2007, 2009; Ranalli 
etal. 2008). We made use of the spatial resolving power of 
Chandra to better inform our estimates of the diffuse emission 
component. To this end, we extracted Chandra “diffuse com¬ 
ponent” spectra from the galaxies after removing the bright 
point-source populations. In this exercise, we excluded events 
from regions that were within 1.5 x the radius of the 90% 
encircled-energy fraction PSFs of all bright Chandra detected 
point sources (see § 4.2 below for details on the Chandra 
detected point-source population) and extracted events from 
within the same apertures used for the galaxy-wide spectra 
discussed above. The resulting diffuse component spectra are 
expected to be dominated by hot gas and unresolved, low lu¬ 
minosity X-ray point-sources. We note that the bright X-ray 
point-source densities are highest in the central regions of the 
galaxies where star-formation is highest. It is in these same re¬ 
gions where we expect to have the hottest X-ray-emitting gas 
associated with the starbursts. For NGC 3256 and NGC 3310, 
the removal of these regions will therefore likely result in 
the removal of a significant fraction of the emission from the 
hottest gas component. 

As a first step, we fit the diffuse, point-source-subtracted 
Chandra spectra using two thermal components (apec in 
XSPEC; Smith etal. 2001) plus a single power-law (account¬ 
ing for unresolved point sources), to determine the tempera¬ 
tures of the galaxy-wide hot gas{3 In the top portion of Ta¬ 
ble 1, we list the best-fit temperatures that result from this pro- 

15 We modeled the point-source-subtracted Chandra spectra using XSPEC 
model tbabs* (apec + tbabs*apec + tbabs*powerlaw) and 
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cess. Similar to the results presented by Jenkins et al. (2004), 
we find best-fitting temperatures of kT ss 0.3 and 0.9 keV for 
NGC 3256; however, we find that much cooler temperatures 
of kT ps 0.2 and 0.3 keV were required for NGC 3310, with 
the latter component being more absorbed and intrinsically 
powerful than the former component. 

Next, for each galaxy, we modeled the full galaxy-wide 
Chandra plus NuSTAR spectrum using the two thermal com¬ 
ponents (with temperatures fixed at the previous best-fit val¬ 
ues) plus a third, hotter (kT £ 1 keV) and more obscured 
thermal component and a non-thermal component to model 
emission from the point source population. For the point 
source spectral component, we tested both a single and broken 
power-law model to see if breaks (e.g., at Em 3-8 keV) were 
preferred, as is most often the case for individual ULX spec¬ 
tra (e.g., Gladstone etal. 2009; Bachetti etal. 2013; Lehmer 
etal. 2013; Rana etal. 2014; Wik etal. 2014a). Since the 
point-source spectra of our galaxies contain the conglomer¬ 
ated summed emission from primarily ULXs, we do not re¬ 
port fits using more complex models that contain detailed 
physical treatment of individual accretion disks and coro- 
nae (e.g., diskbb, diskpbb, diskpn, comptt, eqpair, 
and dkbbf th). In Table 1, we provide the best-fitting param¬ 
eters for fits to the Chandra plus NuSTAR spectra that include 
the two temperature plasma (with fixed temperatures) com¬ 
ponent, plus an absorbed hotter component (in the case of 
NGC 3256 only), plus an absorbed single or broken power- 
law component. In the bottom panels of Figure 4 we display 
the residuals to these fits. 

For NGC 3256, we find that the sum of the diffuse com¬ 
ponent plus a single absorbed (Nu ~ 3 x 10 22 cm -2 ) power- 
law (T ps 2.6) and an absorbed hot (ssl keV) plasma com¬ 
ponent were sufficient to obtain a good fit to the 0.3-30 keV 
spectrum. The absorbed hot plasma component was imple¬ 
mented to fit a notable Fe-line feature in the NuSTAR spec¬ 
trum at «6.4-6.8 keV. Such a feature is also notable in the 
Chandra spectrum of the full galaxy but not the diffuse re¬ 
gion with point sources removed (see above). Further anal¬ 
ysis indicates that the Fe emission photons are concentrated 
in a region within r; 20 arcsec of the nucleus, where most of 
the X-ray detected point sources are located. This component 
was not in our diffuse, point-source-subtracted Chandra spec¬ 
trum because of the spatial coincidence with point sources. 
We note that X-ray binaries and obscured AGN may also pro¬ 
vide some contribution to such a feature. However, given the 
steepness of the galaxy-wide spectral slope above 8 keV, it 
is unlikely that the Fe-line is powered by an obscured AGN, 
which would have a harder spectral slope (F ps 1 . 5-2.0; see 
§4.3). Replacing the single power-law with a broken power- 
law does not improve the quality or change the character of 
the fit (see lower panels of Fig. 4 left). Since the thermal 
components in the galaxy dominates the «0.3-2.5 keV emis¬ 
sion in NGC 3256, a true break in the power-law component 
could be masked if it resides at E 3 keV. 

For NGC 3310, we find that the diffuse component plus a 
broken power-law model provides the best fit to the data; no 
additional absorbed hotter component was required to fit the 
spectrum, as we do not see clear evidence for an Fe-line. The 
best-fitting slopes (Ti ps 2.1 and T 2 ps 2.7) and energy break 
(£Weak ~ 4.7 keV) are within the range of values found for 

fixed the initial absoiption to Galactic and the power-law spectral slope to 
r = 1.8, a value consistent with that of bright X-ray binaries (e.g., Mineo 
etal. 2012a). The apec component abundances were fixed to solar. 


ULXs (e.g., Gladstone et al. 2009). In contrast to NGC 3256, 
which had strong emission features related to hot gas out to 
E ps 3 keV, we find that the overall spectrum of NGC 3310 
above psI keV is dominated by the point-source population. 
As we will discuss in more detail below, the SFR-normalized 
point-source emission in NGC 3310 is a factor of rs 3-10 
times higher compared to other starburst galaxies studied in 
our program (i.e., NGC 253, NGC 3256, and M83). The 
relatively strong emission from the point-source population 
may therefore be masking lower-intensity line features that 
would betray the presence of a hot X-ray plasma or an ob¬ 
scured AGN (e.g., Fe emission-lines). If NGC 3310 harbors 
an obscured AGN, it is almost certainly of low luminosity and 
is likely to provide only a minor perturbation on the non-AGN 
emission in the NuSTAR band. We return to the discussion of 
the potential for AGN in both NGC 3256 and NGC 3310 in 
§4.3 below when we discuss the Chandra properties of the 
nuclear sources in detail. 

4.2. Chandra Point-Source Analyses 

Our joint NuSTAR and Chandra spectral analyses indicate 
that the galaxy-wide 0.3-30 keV emission from NGC 3256 
and NGC 3310 are dominated by the combination of hot gas 
and ULX populations, with no obvious signatures of lumi¬ 
nous AGN. In this section, we gain further insight into the 
X-ray emitting populations in these galaxies by performing 
basic Chandra spectral analyses of the detected point-sources 
and diffuse emission. 

For each galaxy, we identified point sources that had >20 
0.5-7 keV counts that were also detected in the 2-7 keV 
band images using wavdetect at a false-positive probabil¬ 
ity threshold of 10 -5 ; these sources were candidates for per¬ 
forming basic spectral fits. For NGC 3256 and NGC 3310, we 
found 7 and 9 sources, respectively, that satisfied these crite¬ 
ria; the locations of these sources are annotated in Figure 3. 
Taken together, these sources respectively provide «42% and 
rs 79% of the 2-7 keV Chandra net counts within the total 
galaxy apertures defined in §4.1. For these sources, we ex¬ 
tracted the Chandra spectra and fit them using both an ab¬ 
sorbed power-law model (powerlaw) and an absorbed mul¬ 
ticolor accretion disk (diskbb; Shakura & Sunyaev 1973; 
Mitsuda etal. 1984). In this procedure, point-sources were 
extracted from circular regions that encompassed 90% of 
the encircled energy of the PSF (ps2 arcsec in radius for all 
sources). We utilized the same background regions used in 
§4.1 to estimate background spectra for the sources. As such, 
some source regions will contain background related to dif¬ 
fuse emission. Since the exposures are relatively shallow and 
diffuse emission gradients are high, we chose not to attempt 
to model the local background component. When relevant, 
we have qualified interpretations of the Chandra source spec¬ 
tra with this limitation. Given the small number of counts per 
source (33-329 0.5-7 keV counts), all spectral fits were per¬ 
formed by minimizing the Cash statistic (cstat; Cash 1979) 
using spectra binned to a minimum of 1 count per channel. 

In Table 2, we present the best-fit parameters for the point- 
source spectral fits. We note that the parameters in Table 2 
provide only a basic description of the data and are not well 
constrained, e.g., due to degeneracies between column density 
and spectral slope or inner disk temperature. In particular, for 
sources that were located near the nuclei of the galaxies, some 
contributions from hot gas are expected, which will steepen 
the inferred spectral slopes. In the majority of cases, a simple 
power-law provided an acceptable fit to the spectra; we find 
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TABLE 2 

Chandra Point-Source Properties 


powerlaw diskbb 


q J2000 <5j 2000 ^H,int Ab.int kRn log/2-10 keV l°g^-2-10keV Net Counts 

Source ID (hr) (deg) (10 22 cm 4 ) T C (10 22 cm -2 ) (keV) C u (log erg cm -2 s -1 ) (log ergs -1 ) (0.5-7 keV) Note 

(1)(2)(3)(4) (5) (6) (7) (8) (9) (10)(11)(12)(13)(14). 

NGC 3256 

1 . 10 27 50.0 -43 54 19.8 0.02^; 1 2 3 4 5 6 7 8 9 ? L4+^ 43 < 0.16 1.9^; 2 45 32 -13.5 39± 43 ± 7 

2 . 10 27 50.6 -43 54 19.8 <0.38" 1.2f° :2 60 <0.20 2.0)Jj 60 55 -13.2 39.9 67 ± 8 

3 . 10 27 50.8 -43 54 15.2 0.37^ g 2.2^ 44 0.08^ 1.2j° : | 46 52 -13.5 39.7 71 ± 8 

4 . 10 27 51.2 -43 54 13.9 0.36*™ 2.8^ 123 0.04*> :i » 0.8*> :i 126 140 -13.1 40.1 329 ± 18 Nucleus 

5 . 1027 51.6 -43 5409.3 0.32*> : j| 2.9*^ 65 0.02* 1 * 1 0.8*] :2 69 88 -13.5 39.6 139 ± 12 

6 . 10 27 51.8 -43 54 13.3 0.55+g : |f 3.6*>| 60 0.18^g : Jf 0.6)* 2 63 68 -13.7 39.4 118 ±11 

7 . 10 27 55.1 -43 54 46.6 <0.71 0.5^ 33 <0.43 > 1.35 33 25 -13.1 40.0 33 ± 6 

NGC 3310 

1 . 10 38 43.3 +53 31 02.0 0.22^ ]J 2.3^ 31 0.0 1*]-'" 0.9^ 2 27 23 -12.8 308 168 ± 13 

2 . 10 38 44.8 +53 30 04.3 0.55^ 1.6* 1 ; 2 12 0.34+g;f^ 1-T^g-g 11 10 -12.7 40.0 152 ± 12 

3 . 10 38 45.8 +53 30 12.2 0.30* 221 0.4*> :2 9 0.58^g : |g > 10 10 8 -12.4 40.3 174 ± 13 Nuclei] 

4 . 10 3846.0 +53 3004.8 0.71^> ; « 1.8*> : ji 9 0.40* 22 * 1.5^ 10 16 -12.6 40.1 131 ± 11 

5 . 10 38 46.5 +53 30 38.3 0M^ 2.0*} 2 14 0.51:;;' 1.3*” 12 12 -12.7 40.0 142 ± 12 

6 . 10 38 46.6 +53 30 13.7 0.59*>;g 1.7*];jj 2 0.24*>;*} 1.9*$ 3 8 -12.8 39.9 130 ±11 

7 . 103846.7 +53 3007.8 <0.22 1.9*” 5 < 0.08 0.8*” 8 7 -13.1 39.5 213 ± 15 

8 . 10 38 47.2 +53 30 27.9 0.28^;|| 1.8*” 13 0.04*]; 22 1.5*” 13 10 -12.9 39.8 121 ±11 

9 . 10 38 50.2 +53 29 26.0 046^ 2.0/” 22 0.25$$ L0$j 19 12-12.9_308_294 ± 17 


NOTE.—All quoted errors on fit parameters indicate 90% confidence intervals, while quoted errors associated with net counts are 1 a. Col.(l): Source ID. 
Col.(2) and (3): Right ascention and declination, respectively, based on Chandra source location. Col.(4)-(6) provide parameters to spectral fits for a power-law 
with intrinsic absorption model. Col.(4): Intrinsic absorption column density in units of 10 22 cm -2 . Col.(5): power-law photon index. Col.(6): Minimum 
cstatistic C. Col.(7)-(9) provide parameters to spectral fits for a multicolor accretion disk (diskbb) with intrinsic absorption model. Col.(7): Intrinsic 
absorption column density in units of 10 22 cm -2 . Col.(8): Best-fitting inner accretion disk temperature (kT m ) in units of keV. Col.(9): Minimum cstatistic 
C. Col.(10) Number of degrees of freedom for both power-law and multicolor accretion disk fits. Col.(l 1) and (12): Observed 2-10 keV flux and luminosity 
based on the best-fitting model. Col.(13) Net counts in the 0.5-7 keV Chandra band. Col.(14) Notes on the individual sources. 


a slight (~ 1 tr level) statistical preference for the diskbb 
model for sources 1 and 9 of NGC 3310. The luminosity 
range of the point-source populations cover L 2-10 keV = (3- 
20) x 10 39 erg s -1 , indicating that if all sources are discrete ob¬ 
jects, then they are all in the ULX range of luminosities. The 
photon indices for the point sources range from T = 0.4-3.6 
(median (T) = 1.8), with only source 3 in NGC 3310 hav¬ 
ing T < 1.5 at >90% confidence. As we will discuss in the 
next section, source 3 is coincident with the nuclear region in 
NGC 3310 and is a candidate obscured AGN. 

4.3. Contributions from Nuclear Sources and Potential AGN 

Both NGC 3256 and NGC 3310 contain Chandra de¬ 
tected point-like sources coincident with their galactic nuclei 
(NGC 3256 source 4 and NGC 3310 source 3); however, for 
NGC 3256, we only detect the northern nucleus and not the 
southern nucleus. In this section, we investigate the nature of 
the sources coincident with the nuclei using Chandra spec¬ 
tral constraints and discuss their potential contributions to the 
0.3-30 keV spectra. 

For NGC 3256, there is a clear detection of the northern nu¬ 
cleus, which is the brightest point-like source in our Chandra 
catalog (source 4). Using a deeper «30 ks Chandra exposure 
(ObsID 835), Lira etal. (2002) found that the nuclear emis¬ 
sion is clearly extended, suggesting that the emission must 
include extended and/or multiple point-sources. From Ta¬ 
ble 2, we see the Chandra spectrum of this source is fit well 
with a steep power-law (T = 2.8 ± 0.4) with modest absorp¬ 
tion (Ah ps 4 x 10 21 cm -2 ). Such a steep photon index is un- 


characteristc of AGN, which typically span T p> 1-2.2. We 
refit the spectrum using only the 3-7 keV range, to mitigate 
any biases in the spectral fit due to soft emission from the 
hot gas, which will inevitably be coincident with our spectral 
extraction of this region. The resulting best-fit photon index 
remained unchanged, indicating that there is no evidence for 
an emerging hard X-ray AGN in this region. This result is 
consistent with previous studies that have concluded that the 
northern nucleus in NGC 3256 is powered almost exclusively 
by star-formation activity (e.g.. Lira etal. 2002; Jenkins etal. 
2004; Sakamoto etal. 2014). In Figure 4, we show the con¬ 
tribution from the northern nucleus extrapolated through to 
the NuSTAR band. Based on this extrapolation and the global 
fit presented in §4.1, we estimate that this source will con¬ 
tribute f>;4—10% of the galaxy-wide emission across the full 
0.3-30 keV band. 

Our 16 ks Chandra observation of NGC 3256 did not yield 
a clear detection of the southern nucleus (location indicated 
in the left panel of Fig. 3), which has previously been shown 
to have submm and near-infrared evidence for a collimated 
outflow, potentially due to an AGN that is currently inactive 
(e.g., Sakamoto etal. 2014; Emonts etal. 2014). Two pre¬ 
vious deeper «25 ks archival Chandra exposures (ObsIDs 
835 and 3569) reveal clear detections of a source coincident 
with the southern nucleus, suggesting that the source does 
produce powerful, albeit heavily obscured. X-ray emission 
at the «10 40 erg s -1 level (in the 0.3-10 keV band). Us¬ 
ing the Chandra exposure from ObsID 835, as well as radio, 
near-IR, and optical data. Lira etal. (2002) showed that this 
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source is likely to be a low luminosity AGN that does not 
provide a significant contribution to the total X-ray emission 
of the galaxy. No evidence was found implicating a lumi¬ 
nous, heavily obscured AGN (e.g., Fe-line emission). The 
non-detection of the southern nucleus in our 16 ks Chandra 
exposure indicates a change in the southern source with re¬ 
spect to the previous two observations. In 2000 Jan 5 and 
2003 May 23, the southern source was detected with similar 
2-7 keV Chandra count-rates of (2.7 ±0.3) x 10 -3 cnts s -1 
and (2.9 ± 0.3) x 10 -3 cnts s -1 , respectively. Whereas the 
3cr upper limit on the 2-7 keV Chandra count-rate in our 
2014 Aug 28 observation is <1.2 x 10~ 3 cnts s -1 , indicating 
a factor of £ 2.0-2.5 times decline in the southern source 
hard-band intensity. As shown by our NuSTAR spectrum, 
the steep spectral slope of the full galaxy persists out to at 
least «30 keV, limiting the possible contributions from a hard 
southern nuclear AGN. This implies that if there is an AGN 
associated with the southern nucleus, it is likely to be of low 
luminosity. 

To place more stringent constraints on the spectral prop¬ 
erties of the southern nucleus, we combined the previous 
two Chandra observations of NGC 3256 to create a merged 
«50 ks data set. We extracted the merged spectrum of 
the southern nuclear source and performed spectral fitting in 
XSPEC. The spectrum can be fit well by a simple power- 
law model with an inverted photon index of T rs —0.1; the 
best-fit model implies a 2-10 keV luminosity of L 2-10 kev ~ 
10 4() erg s -1 . It is likely that the spectrum contains both contri¬ 
butions from hot gas and a heavily obscured power-law source 
coincident with the nuclear region. The black hole mass 
for the southern nucleus is on the order of ~ 10 7 M & (e.g., 
Alonso-Herrero et al. 2013), implying an Eddington luminos¬ 
ity of ~10 45 erg s -1 . If the southern nucleus is indeed an AGN 
it would therefore have T 2-10 kev/CEdd ~ 10 -5 , and the ex¬ 
pected intrinsic spectral slope of this source would be between 
T = 1.5-2.0 (Younes etal. 2011). Using these constraints, we 
fit the spectrum of the southern nucleus using an absorbed 
plasma plus an absorbed power-law model (tbabs*apec+ 
tbabs*powerlaw in XSPEC). We held the plasma temper¬ 
ature fixed at 0.6 keV (see Lira et al. 2002 for motivation) and 
the power-law slope fixed at T = 1.8 and fit for the component 
normalizations and obscuring columns. This model yielded a 
good fit to the data for column densities iV H gas rs 10 22 cm -2 
and Nh.pl ~ 5 x 10 22 cm -2 . 

In the left panel of Figure 4, we show the Chandra spectrum 
and best-fitting model of the southern nucleus of NGC 3256 
from the previous two observations when the source was de¬ 
tectable ( orange data and cun’e). The data and model have 
been scaled down by a factor of 2.5 to show the upper limit 
to our 2014 observation (see above). Based on an extrapola¬ 
tion of this model to the NuSTAR band, the southern nucleus 
is expected to contribute <, 10-30% of the NuSTAR emission 
at 10-30 keV. Such a low level of AGN activity is insuffi¬ 
cient to have any material effect on the overall shape of the 
galaxy spectrum nor on the constraints presented in §4.1 and 
Table 1. We therefore expect that the 0.3-30 keV spectrum of 
NGC 3256 is dominated by non-AGN emitting sources (i.e., 
hot gas and X-ray binaries). 

For NGC 3310, a previous «50 ks Chandra observation 
from 2003 (ObsID 2939) was used by Tzanavaris & Georgan- 
topoulos (2007) to show that there is evidence for an Fe-Ka 
emission line in the nuclear source (source 3), implying the 
presence of a hidden AGN. Our r; 10 ks Chandra observation 
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Fig. 5.— NuSTAR light profiles for NGC 3310 in the 4-8 keV (top) and 
8—25 keV ( bottom) bands. A blurred PSF model (blue dashed curves) and 
a model based on the relative 4-8 keV Chandra point-source brightnesses 
(“Chandra Model”; red solid curves) have been compared with the NuSTAR 
light profiles (black points). We find that the “Chandra Model” characterizes 
well the NuSTAR galaxy light profile in the Chandra and NuSTAR overlap¬ 
ping 4-8 keV bandpass. However, due to the presence of harder sources in 
the central region of the galaxy, as compared with sources in the outskirts, 
the 8-25 keV light profile becomes more point-like. Such a result could be 
explained by the presence of an obscured AGN and/or harder ULXs in the 
central region of the galaxy. 

was too shallow to verify the claimed Fe-Ka feature, which 
was of low significance («2 cr) in the much deeper 2003 ob¬ 
servation. However, we do find that the best-fitting power-law 
spectral slope of the nuclear point-source T « 0.4 is much 
shallower than the other sources in the galaxy (F « 1.6-2.3), 
and the luminosity of the source (L 2-10 kev ~ 10 4() erg s -1 ) is 
the highest of all the point sources. These facts are consis¬ 
tent with an obscured or low luminosity AGN hypothesis. We 
note that if we invoke a more complex model appropriate for 
AGN that includes both direct and scattered components (e.g., 
pexrav or MyTorus; Magdziarz & Zdziarski 1995; Mur¬ 
phy & Yaqoob 2009), the expected contribution to the galaxy¬ 
wide emission above 10 keV can quickly approach 100%, and 
the intrinsic emission could be much higher than the observed 
luminosity inferred by Chandra. 

Following a similar approach to that taken for the southern 
nucleus in NGC 3256, we tested the influence that the nuclear 
source in NGC 3310 would have on the ss7-30 keV emission, 
if it were truly an AGN, by extrapolating the Chandra spec¬ 
trum through to the NuSTAR band. The black hole mass of 
NGC 3310 is currently poorly constrained, but has been esti¬ 
mated via gas kinematics to be ~(5—40) x 10 6 M 0 , which im¬ 
plies an Eddington luminosity of ~(6-50) x 10 45 erg s -1 (Pas- 
torini etal. 2007). The 2-10 keV luminosity of the source, as 
measured by Chandra, is L 2 -iokeV ~ 10 4l) erg s -1 , implying 
T 2-10 kev/CEdd ~ 10 -6 — 10 —5 . From Younes etal. (2011), this 
implies an intrinsic photon index of T = 1.5-2.2. Following 
the approach used for the southern nucleus of NGC 3256, we 
fit the Chandra spectrum of the NGC 3310 nuclear source us¬ 
ing an absorbed plasma plus an absorbed power-law model 
(tbabs*apec+ tbabs*powerlaw in XSPEC) with the 
plasma temperature fixed at 0.3 keV and the power-law slope 
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fixed at r = 1.8. These parameters yielded a good fit to 
the data for column densities AH,gas ~5x 10 21 cm -2 and 
Ah.pl « 3 x 10 22 cm -2 . 

In the right panel of Figure 4, we show the Chandra spec¬ 
trum for the nuclear source of NGC 3310 and the best-fit 
model extrapolated across the full 0.3-30 keV band. It is 
clear that the nuclear source is expected to provide a non- 
negligible contribution to the X-ray emission above 4-5 keV. 
Based on our extrapolated model, the point-source makes a 
fractional contribution of 3% of the galaxy-wide emis¬ 
sion at E < 2 keV and reaches a maximum contribution of 
~30-100% in the 20-30 keV band. If we revist our fits to the 
galaxy-wide spectrum performed in §4.1 and now fix the dif¬ 
fuse component and best-ftting model for the nuclear source, 
we find that the remaining spectrum associated with the bright 
point sources can be well characterized by a broken power- 
law with parameters consistent with those presented in Ta¬ 
ble 1. 

4.4. Spatially Extended NuSTAR Emission for NGC 3310 

Although the NuSTAR PSF is too large to identify individual 
X-ray point sources in the galaxies studied here, it is sufficient 
for identifying extended features due to the «50 arcsec wide 
distribution of point sources throughout NGC 3310. To test 
for extended NuSTAR emission, we first constructed a spa¬ 
tial model of the expected NuSTAR emission profile within a 
bandpass that overlapped with Chandra (i.e., 4-8 keV), and 
then tested to see how the NuSTAR emission profile changed 
in a bandpass above the Chandra response (i.e., 8-25 keV). If 
the X-ray SEDs of the underlying point sources were all sim¬ 
ilar, we would expect the profile to remain the same for the 
soft and hard bands. Using the positions and relative 4-8 keV 
Chandra count-rates for the point sources, we constructed a 
model distribution map of NuSTAR counts. This map was 
constructed by co-adding NuSTAR PSF images centered at the 
locations of the Chandra point sources, with normalizations 
scaled by the 4-8 keV Chandra count-rates. Throughout this 
process, our NuSTAR PSFs were constructed by accounting 
for the off-axis angles of the sources and include errors asso¬ 
ciated with variations in off-axis angle due to the motions of 
the mast. We note that some relative astrometric offset is ex¬ 
pected to be present between the NuSTAR and Chandra data. 
In order to estimate this error, we fit Gaussian profiles to both 
the NuSTAR 4-8 keV image emission and the Chandra- based 
model distribution map to identify the location of the peak 
emission. We then extracted the emission profiles for the im¬ 
ages and models and compared them. 

In Figure 5, we show the spatial distribution of NuSTAR 
counts for the 4-8 keV and 8-25 keV bands. The expected 
distributions for the Chandra- based model (described above) 
and an individual point-source PSF model have been overplot¬ 
ted for comparison. By construction, the 4-8 keV NuSTAR 
emission profile follows the Chandra- based model well and 
is statistically inconsistent with a single point source. How¬ 
ever, at 8-25 keV, the NuSTAR emission profile is more con¬ 
sistent with a single point-source (i.e., the NuSTAR PSF) than 
with the Chandra- based model, implying that the central re¬ 
gion of NGC 3310 contains a harder X-ray spectrum (across 
8-25 keV) compared to point-sources in the outskirts of the 
galaxy. Our spectral extrapolation of the nuclear source pre¬ 
sented in §4.3 and Figure 4 (right) shows that the fractional 
contribution that the source makes to the galaxy-wide emis¬ 
sion increases with energy. It is also the case that the two 


sources in the outskirts of the galaxy (sources 1 and 9; see 
Fig. 3 and Table 2) have the steepest spectral slopes out of 
all nine point sources. Therefore, the hard spectra of the nu¬ 
clear source and the five relatively hard non-nuclear sources 
that are within r; 20 arcsec of the galactic center are likely to 
be responsible for the transition to a more point-like emission 
profile with increasing energy. As we discussed above, we 
find a steep T « 2.7 power-law slope at E A 4 keV for the 
galaxy-wide spectrum of NGC 3310, which is steeper than 
the best fit single power-law slopes (T rj 0.3-2.0) derived for 
all five sources. This implies that, although the spectra of 
these sources may be intrinsically hard in the Chandra band, 
spectral turnovers at E 3-10 keV must be typical of these 
sources, making it difficult to quantify the relative influences 
these sources have on the emission measured in the NuSTAR 
band. 

5. DISCUSSION 

We have utilized nearly simultaneous Chandra and NuS¬ 
TAR observations of the powerful local starburst galaxies 
NGC 3256 and NGC 3310 to investigate the nature of 
the X-ray-emitting components across the broad 0.3-30 keV 
bandpass. Similar to previous XMM-Newton and Chandra 
studies, we find significant kT r* 0.2-1 keV plasma emission 
that dominates the spectra of these galaxies out to «3 keV 
and ssl keV for NGC 3256 and NGC 3310, respectively. At 
the higher energies probed by NuSTAR (E If 3 keV), we find 
the majority of the emission is produced by populations of 
the brightest 5-10 ULXs. The cumulative spectra of both 
galaxies exhibit steep power-law slopes (T « 2.6-2.7) above 
~3-6 keV that are similar to those found in NuSTAR stud¬ 
ies of individual UFXs, implying that no strong hard com¬ 
ponent exists for either galaxy. This places stringent limits 
on the influence of possible AGN candidates in both galax¬ 
ies, which are expected to have much flatter spectral slopes 
in this band. Sensible extrapolations of the Chandra spectra 
of AGN candidates in these galaxies, the southern nucleus of 
NGC 3256 and the nucleus of NGC 3310, indicate that these 
candidates provide only minority contributions to the galaxy¬ 
wide NuSTAR emission out to 30 keV. Regardless of extrap¬ 
olation assumptions, our constraints indicate that the super- 
massive black holes in these galaxies are of low intrinsic lu¬ 
minosity and are accreting at the 10~ 6 -10~ 5 Eddington lev¬ 
els. Our spectral constraints on NGC 3256 and NGC 3310 
therefore provide first measurements of the 0.3-30 keV spec¬ 
tra of the non-AGN X-ray emitting populations within power¬ 
ful starburst galaxies with high-sSFRs. 

In Figure 6, we show the SFR-normalized best-fitting 
0.3-30 keV SEDs (i.e., uLJSFR) of NGC 3256 and 
NGC 3310. For comparison, we also provide best-fit mod¬ 
els for M83 (Yukita etal. in-preparation) and NGC 253 (Wik 
etal. 2014a), which were also constrained using nearly simul¬ 
taneous observations with Chandra/XMM-Newton and NuS¬ 
TAR. The displayed SEDs have not been corrected for intrin¬ 
sic absorption and are therefore representative of the observed 
X-ray spectra. In the low-energy regime (E f 2-3 keV), all 
four galaxies have significant contributions from line-emitting 
hot plasmas. With the exception of NGC 253 ( orange curve) 
the SFR-normalized gas components have similar normaliza¬ 
tions around r ;1 keV, consistent with previous studies that 
have found a direct scaling of the hot gas emission with SFR 
(see, e.g., Mineo et al. 2012b). Compared with the other three 
galaxies, NGC 253 is relatively edge-on (i Ri 80 deg), and 
therefore absorption is likely playing a role in the apparent 
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Fig. 6.— SFR-normalized 0.3-30 keV SEDs for NGC 3256 (red), NGC 3310 (green), NGC 253 (orange), and M83 (blue). For illustrative purposes, we 
have displayed only the best-fit SEDs to the nearly simultaneous Chandra/XMM-Newton plus NuSTAR data. NGC 253 and M83 fits are taken from the best-fit 
models presented in Wik etal. (2014a) and Yukita etal. (in-preparation). For comparison, we have shown NuSTAR constrained SEDs for the ULX NGC 1313 
X-l (dotted cun’e; Bachetti etal. 2013), the luminous black-hole X-ray binary GRS 1915+105 in the low/hard state (short-dashed curve'. Miller etal. 2013), and 
the Compton-thick AGN in Arp 299 (long-dashed curve; Ptak etal. 2014); these SEDs have been normalized to vL ^/SFR = 1.3 X 10 39 erg s 1 (Mq yr 1 )" 1 at 
5 keV for ease of comparison. In this representation, a power-law SED with photon-index T = 2 is a flat line. It is clear that all four galaxies here show steepening 
in their SEDs above Ri3-6 keV, where their spectral slopes reach T ~ 2.5-3.0 due primarily to the ULX population. These constraints indicate that the ULX 
population on average has a similar spectral shape to that of super-Eddington accreting sources like NGC 1313 X-l. 


deficit of hot gas emission in the galaxy. 

At higher energies, E k, 3-30 keV, we find that all galax¬ 
ies, except for NGC 3310 (see discussion below), have SFR- 
normalized SEDs that are in good agreement with each other. 
For all four galaxies, the steepening of the spectral slope in the 
~6-30 keV range is clearly apparent. For M83, it appears that 
the spectral slope in this regime is somewhat steeper than it is 
for the other three galaxies. As discussed by Yukita etal. (in¬ 
preparation), the spectral shape of M83 is strongly influenced 
by a single variable ULX (see also Soria et al. 2012 for details 
on this source). For comparison, we have overlaid in Figure 6 
the NuSTA R-constrai ncd SEDs (scaled to the mean starburst 
galaxy emission at 5 keV) of the ULX NGC 1313 X-l (Ba¬ 
chetti et al. 2013), the luminous black-hole X-ray binary GRS 
1915+105 in the low/hard state (Miller etal. 2013), and the 
Compton-thick AGN in Arp 299 (Ptak etal. 2014). From this 
illustration, it is clear that the > 10 keV starburst galaxy SEDs, 
including the turnover in spectral slope, are very similar to 
the SED of NGC 1313 X-l, a likely super-Eddington accret¬ 
ing ULX. Given the fact that ULXs clearly dominate the point 
source emission from this galaxy population, our observations 
constrain the average spectral shape for the ULX population, 
in general. 

Interestingly, for NGC 3310, we observe a clear factor of 
rG- 10 times excess X-ray emission per unit SFR over the 
~6-30 keV range. The excess is the result of an overabun¬ 
dance of ULXs in the galaxy compared to typical galaxies. 
FromMineo etal. (2012a), we expect that there would be sal- 
4 ULXs given the SFR of NGC 3310 (SFR « 6 M 0 yr 1 ); 
however, 9 ULXs are clearly detected, indicated a significant 


excess in both the X-ray emission and the number of ULXs. 
In general, an excess of luminous X-ray binaries can be ex¬ 
plained by either (1) a star-formation history that is heavily 
weighted towards an epoch when ULXs are expected to be 
most luminous due to the presence of the most massive donor 
stars («5-10 Myr) and/or (2) star-formation that is happening 
in low-metallicity environments, in which an excess popula¬ 
tion of black holes form (e.g., Linden et al. 2010). 

//5T-based studies of the young star cluster properties in 
NGC 3310, using SED fitting of UV-to-near-IR HST data, 
indicate a peak intensity in the recent star-formation his¬ 
tory around rGO Myr ago (de Grijs etal. 2003a,b), well 
past the peak of ULX activity. Furthermore, similar star- 
cluster ages are also seen in NGC 3256 (e.g., Trancho et al. 
2007), which suggests that the excess of luminous X-ray 
binaries in NGC 3310 is unlikely to be caused by differ¬ 
ences in recent star-formation history alone. By contrast, the 
metallicity distribution of young star clusters in NGC 3310 
has been estimated to peak at Z k 0.4 Zq (de Grijs etal. 
2003a,b), compared with Z w 1.5 Z 0 for M83, NGC 253, 
and NGC 3256 (e.g., Zaritaas etal. 1994; Boselli etal. 2002; 
Bresolin etal. 2014; Trancho etal. 2007). The relatively 
low metallicity in NGC 3310 is the result of the recent 
cannibalization of a low-metallicity dwarf galaxy that trig¬ 
gered the current star-formation event. It is therefore plau¬ 
sible that the lower metallicity of NGC 3310 is the driv¬ 
ing mechanism behind the excess of ULXs per unit SFR 
in this galaxy. Indeed, past studies of both the number 
of ULXs and galaxy-wide X-ray luminosity per unit SFR 
have revealed excesses at low metallicities (e.g., Mapelli 
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etal. 2009; Prestwich etal. 2013; Basu-Zych etal. 2013b; 
Brorby etal. 2014). From the population synthesis models 
from Fragos etal. (2013b), galaxies with Z ps 1.5 Z 0 and 
0.4 Z 0 are predicted to have SFR-normalized 2-10 keV lu¬ 
minosities of Li- iokev/SFR = (1.1-3.4) xlO 39 erg s” 1 and 
(3.5-11.2) xlO 39 erg s -1 , respectively (e.g., see their Fig. 2). 
We measure corresponding values of ^ 2-10 kev/SFR = 2.1 x 
10 39 erg s _1 and 7.3 x 10 39 erg s -1 , respectively, in very good 
agreement with the theoretical model predictions. We there¬ 
fore favor metallicity as being the dominant underlying factor 
responsible for the observed p» 3—10 times excess X-ray emis¬ 
sion per unit SFR in NGC 3310; however, future investiga¬ 
tions of the properties of the immediate environments (i.e., 
characteristic stellar ages and metallicities) near the ULXs 
would help to better discriminate the roles that age and metal¬ 
licity play in producing the excess ULX population. 
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